Report  1539 


m  -i  zp  ^ 

O  i  \{  pc 


USAOAC  TECHNICAL  LIBRARY 


:  J;I’!*,*’v’*,*V** 


?.*.•.  : :  ::V:  v.v  •*.*.* 


5  0712  01017315  0 


r.v  *;*.*.• 


_ _ _ _ _ _ _ _ 


AERODYNAMICS 


STRUCTURAL 

MECHANICS 


APPLIED 

MATHEMATICS 


I 


c  1 


•  \  V 


HYDROMECHANICS  I  ^ 


BLAST  TRANSMISSION  INTO  CHAMBERS 


:nj? 

Sft 

ift 


LIBRARY  CORY 

return  to 

OCT  2  5  W4 


by 


George  Chertockf Ph.D. 


NAVAL  UKLlif\r»Ct 

Ilaeoratory  library 


LIBRARY  COPY 
RETURN  TO 

OCT  2  9  1974 

NAVAL  OrvuwntE 
LABORATORY  LIBRARY 


STRUCTURAL  MECHANICS  LABORATORY 

RESEARCH  AND  DEVELOPMENT  REPORT 


November  1961 


Report  1539 


PRMC-THB-6*8  (Rev.  3-58) 


BLAST  TRANSMISSION  INTO  CHAMBERS 


by 


George  Chertock,Ph.D. 


November  1961 


Report  1539 


TABLE  OF  CONTENTS 


Page 


ABSTRACT  .  1 

INTRODUCTION  .  1 

FLOW  EQUATIONS  .  1 

STEADY  FLOW  RESISTANCE  .  8 

RANGE  OF  APPLICATION  .  10 

APPENDIX  -  SAMPLE  CALCULATION  .  12 

REFERENCES  .  14 


LIST  OF  FIGURES 


Figure  1 
Figure  2 

Figure  3 

Figure  4 


Page 

Variation  of  Chamber  Pressure  with  Time 

for  =  2.0,  t=  0.5,  and  k  =  0,  0.12,  and  0.20...  5 

Dependence  of  Maximum  Overpressure  in  Chamber 


on  Fill  Time  of  System  for  k  =  0  and  Blast 

Pressures  of  0.25,  0.50,  1.0,  2.0,  and 

4.0  Atmospheres  .  6 

Variation  of  Maximum  Overpressure  in  Chamber 

with  k,  t,  and  y^Q .  7 

Duct-Chamber  System .  12 


ii 


ABSTRACT 


Equations  are  derived  for  the  pressure  rise  in  a 
terminal  room  or  chamber  due  to  the  propagation  through 
a  duct  system  of  a  transient  pressure  from  an  external 
blast  pressure  wave.  Numerical  values  for  the  peak 
pressure  in  the  chamber  are  given  in  terms  of  the 
parameters  of  the  blast  pressure  wave  and  the  character¬ 
istics  of  the  duct-chamber  system.  The  calculation 
of  the  relevant  characteristics  of  the  duct-chamber 
system  is  described  and  illustrated  in  an  appendix. 


INTRODUCTION 

A  significant  hazard  from  an  air  burst  of  a  nuclear  weapon  arises 
from  the  possible  penetration  of  the  blast  pressures  into  a  structure. 

The  propagation  of  relatively  low  air-blast  pressures  down  a  smokestack 
could  severely  damage  a  fire-box  or  boiler  system.  The  penetration  of 
blast  pressures  through  a  ventilation  duct  and  into  a  protective  shelter 
could  do  serious  damage  within  the  shelter,  yet  the  same  pressures  might 
be  ineffective  against  the  outside  of  the  shelter.  In  order  to  devise 
protective  measures  against  such  hazards,  it  is  desirable  to  analyze 
the  propagation  of  these  blast  pressures  and,  in  particular,  to  estimate 
the  peak  pressure  transmitted  to  an  interior  chamber. 

The  anlaysis  is  based  on  a  plausible  generalization  of  the  equations 
of  steady  flow  through  the  duct  and  uses  gross  parameters  of  the  duct- 
chamber  system  which  can  easily  be  evaluated.  It  is  therefore  applicable 
to  practical  situations  involving  real  structures.  An  exact  analysis 
might  conceivably  be  made  by  a  solution  of  the  exact  flow  equations  for 
every  point  on  the  duct  at  every  instant,  but  it  seems  evident  that 
this  method  would  be  prohibitively  difficult  except  perhaps  for 
exceedingly  simple  geometries. 

FLOW  EQUATIONS 

In  this  analysis  the  distributed  flow  resistance  in  the  duct  system 
is  divided  into  two  components;  a  lumped,  steady  flow  resistance,  and 
a  lumped  inertial  resistance.  The  steady  flow  resistance,  is  assumed  to 
be  proportional  to  the  mean  instantaneous  dynamic  pressure  in  the  duct, 
with  the  proportionality  constant  depending  only  on  the  shape  of  the 


duct  system.  The  justification  for  this  assumption  will  be  discussed 
in  a  later  section.  The  inertial  resistance  is  taken  as  equal  to  the 
rate  of  change  of  the  total  momentum  of  air  in  the  duct. 

Thus,  if  the  approach  velocity  of  the  air  outside  the  duct 
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is  negligible  and  if 

p^  is  the  pressure  in  the  external  blast  wave, 

P^>  the  pressure  in  the  terminal  chamber, 

P,  the  average  air  density  in  the  duct, 
v,  the  average  air  velocity  along  the  duct, 

1,  the  length  of  the  duct, 

A,  the  cross-sectional  area  of  the  duct  where  velocity  is  v, 
q  =  pAv,  the  mass  flow  rate  in  the  duct, 
c,  the  "loss  coefficient"  for  the  duct,  and 
V,  the  volume  of  the  terminal  chamber, 

then 


P 


1 


v 


d( Plv) 
dt 


cqjq  1  dq 

—  2  + - 

2  PA  A  dt 


Cl] 


The  mean  density  in  the  duct  is  assumed  to  be  equal  to  the  average 
of  the  densities  at  the  two  ends,  i.e., 


2P  =  Pi  +  P2 


[2  ] 


The  pressure  at  the  entrance  to  the  duct  is  assumed  to  vary  with  time 
in  the  same  way  as  a  standard  free-field  blast  wave,  and  the  density 
is  assumed  to  decay  isentropically  from  its  initial  shock  w'ave  peak. 


2 


thus 


Pj  =  P0  +  (PlO  -  P0)  e't/e  (1-|)  [3] 

P_  _  6  p10  +  PQ  /  P1  \ 

po  P10  +  6  po  v  pio/  t4  3 


where  and  9  are  the  peak  pressure  and  the  duration  of  the  incident 
blast  wave  and  p^  and  PQ  are  atmospheric  pressure  and  density.  The  total 
mass  of  air  in  the  terminal  chamber  must  increase  at  the  rate  of  the 
mass  inflow  into  the  chamber. 


d_ 

dt 


[5] 


The  mean  pressure  in  the  terminal  chamber  is  assumed  to  rise  isentropically . 
This  assumption  is  a  matter  of  computation  convenience  and  neglects  any 
effects  due  to  the  delay  in  bringing  the  inflow  to  rest  and  the  dissipation 
processes.  Hence 


Equations  [1  ]  through  [6  ]  are  sufficient  to  specify  completely  the 
time  variation  of  the  pressure  and  density  in  the  terminal  chamber  as 
a  function  of  six  parameters,  namely,  the  peak  pressure  and  duration  of 
the  blast  wave,  the  shape,  length  and  volume  of  the  duct  and  the  volume 
of  the  terminal  chamber.  However,  by  a  proper  choice  of  nondimens ional 
units,  the  six  parameters  can  be  combined  into  three  nondimens ional 
parameters  which  completely  specify  the  problem  in  a  more  economical  way, 
Thus,  if  a  new  set  of  units  for  pressure,  density,  and  time  is  defined  by 
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y=l 


z  = 


x '  "§■ 


[7] 


the  equations  reduce  to 


2  ^  z2 

t~rL  + 

dx 


4-  za  \  dx  / 


Z1  +  Z2 


dz. 


dx 


=  y0  -_Jl 

(YT)  2 


[8] 


yl  =  1  + 


(>•10  -  1 )  (1-) 


[9] 


ll  “ 


6  y10  +  1 
y10  +  6  V  y10 


5/7 


[10] 


y2  =  Z2 


7/5 


[11] 


where,  there  remain  only  three  nondimens ional  parameters  y^g,  T>  and  k. 
The  quantity  y^g  is  the  peak  pressure  in  the  external  blast  wave,  measured 
in  units  of  atmospheres .  The  parameter  t  is  defined  by 


[12] 


and  may  be  interpreted  as  a  nondimensional  "fill  time"  (fill  time  in 
units  of  the  duration  0  )  for  the  duct-chamber  system.  Note  that  t 
depends  only  on  the  steady  flow  characteristics  of  the  duct  and  is 
independent  of  the  duct  length. 
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The  adiabatic  constant  V  in  Equation  [8  ]  could  have  been  incorporated  into 
the  definition  of  t,  but  it  was  desired  to  retain  the  same  definition 
as  previously  used.1  The  remaining  parameter  k  is  a  measure  of  the 
relative  importance  of  the  terms  for  unsteady  flow  resistance  and  for 
steady  flow  resistance  and  is  defined  by 


=  i /3I 

V  cV 


[13] 


2 

where  A  here  is  the  mean  sectional  area  along  the  duct.  Thus  k  is 
directly  proportional  to  the  ratio  of  duct  volume  to  chamber  volume  and 
ig  inversely  proportional  to  the  loss  coefficient  in  the  duct. 

Equations  [8]  through  [11  ]  are  a  set  of  nonlinear  equations  which 
must  be  solved  numerically  for  each  set  of  parameters.  This  can  easily 
be  done  with  a  high-speed  computer,  and  calculated  values  of  the  peak 
chamber  pressure  have  been  compiled  for  a  wide  range  of  values  of  the 
parameters.  In  fact,  a  program  is  available  for  the  solution  of  the 
problem  on  the  IBM  7090  computer  for  any  set  of  values  of  y1Q,  t,  and  k 
Typical  solutions  for  the  pres-  l0 
sure  in  the  terminal  chamber  as  a 


function  of  time  are  shown  in 
Figure  1.  These  curves  were  cal¬ 
culated  for  the  special  case  of 
y.^  =  2  (peak  blast  overpressure 
=  1  atm),  t  =  0.5,  and  for  three 
different  values  of  k;  k  =  0,  0.12, 
and  0.20.  When  the  unsteady  flow 
resistance  is  negligible  (k  =  0), 
the  flow  into  the  terminal  chamber 
continues  until  the  chamber  pressure 
is  equal  to  the  external  pressure. 
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Figure  1  -  Variation  of  Chamber 
Pressure  with  Time  for 

y10  =  2'°’  T  =  0,S’ 
and  k  =  0,  0.12,  and  0.20 
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At  that  time  the  chamber  pressure  is  a  maximum,  and  the  flow  reverses. 

The  smaller  the  fill  time  of  the  system,  the  sooner  the  terminal  pressure 
reaches  its  maximum,  and  the  closer  this  maximum  pressure  is  to  the  peak 
blast  pressure..  When  the  unsteady  flow  resistance  is  significant  (k>0), 
the  chamber  pressure  rises  at  a  slower  rate  and  continues  to  rise  after 
the  time  when  it  is  equal  to  the  external  blast  pressure  because  of  the 
inertia  of  the  air  in  the  duct.  The  chamber  pressure  may  finally  reach 
a  peak  value  that  is  not  very  different  than  that  with  k  =  0.  If  the 
characteristic  fill  time  t  is  small  enough  and  k  is  large  enough’,  the 
peak  chamber  pressure  can  exceed  the  peak  blast  pressure. 

The  manner  in  which  the  peak  chamber  pressure  varies  with  the  three 
characteristic  parameters  is  shown  in  Figures  2  and  3,  which  have  been 
compiled  from  the  numerical  solutions.  In  many  practical  situations, 
the  inertial  ratio  k  has  only  a  small  effect  on  the  peak  chamber  pressure, 
and  it  is  particulary  significant  to  consider  first  how  the  peak  over¬ 
pressure  in  the  chamber  varies  with  T  and  with  y^Q  for  k  =  0. 

The  variation  of  the  peak  overpressure  in  the  chamber  is  shown  in 
Figure  2.  The  quantity  plotted  is  the  maximum  value  of  (y2  -  1)  /  (y1Q  -  1) 


Figure  2  -  Dependence  of  Maximum  Overpressure  in  Chamber  on  Fill  Time  of  System 
for  k  =  0  and  Blast  Pressures  of  0.25,  0.50,  1,0,  2.0,  and  4.0  Atmospheres 


6 


Peak  Chamber  Overpressure  {  Atmospheres ) 


Figure  3  -  Variation  of  Maximum  Overpressure  in  Chamber 
with  k,  t,  and  y. 
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which  is  the  ratio  of  the  peak  overpressure  in  the  chamber  to  the  peak  over¬ 
pressure  in  the  blast  wave.  This  method  of  plotting  reduces  the  dependence 
of  the  results  on  the  magnitude  of  the  peak  blast  pressure  and  makes  it  easy 
to  interpolate  the  data  for  values  of  the  blast  pressure  other  than  the  five 
values  calculated.  The  calculations  cover  the  range  from  t  =  0  to  5.0,  and 
for  a  peak  blast  overpressure  up  to  4  atm  (60  psi). 

The  effects  of  all  parameters  on  the  peak  chamber  pressure  are 
shown  in  Figure  3,  which  is  a  plot  of  the  peak  overpressure  in  the  terminal 
chamber  against  the  nondimens ion al  fill  time  t,  with  peak  blast  pressure  and 
inertial  constant  as  variable  parameters.  It  xs  clear  from  the  curves  that, 
as  k  increases  from  0  to  about  0.2,  the  peak  chamber  pressure  also  increases. 
However,  for  higher  values  of  the  inertial  ratio,  the  net  effect  of  k  can 
be  to  decrease  the  peak  chamber  pressure.  Because  of  the  nonlinearity  of 
the  flow  equations,  the  peak  overpressure  in  the  terminal  chamber  is  not 
simply  proportional  to  the  peak  blast  pressure.  The  lower  the  peak  blast 
pressure,  the  higher  is  the  ratio  of  peak  chamber  overpressure  to  peak  blast 
overpressure,  and  the  more  this  ratio  is  affected  by  the  unsteady  flow  term. 


STEADY  FLOW  RESISTANCE 

The  numerical  calculations  show  that  the  most  important  parameter 
in  determining  the  peak  chamber  pressure  is  the  steady-flow-resistance 
term,  and  it  is  important  to  justify  the  form  of  this  term  in  the  equations. 
For  this  purpose,  consider  a  steady  flow  between  two  ends  of  a  complex 
duct,  where  the  area  at  the  ends  is  large  enough  so  that  the  dynamic 
pressure  can  be  neglected  in  comparison  with  the  dynamic  pressure  in  the 
duct.  The  flow  equation  becomes 

2 

[14] 


> 

c  q" 


Pl  -  P2  = 


2  P  A 


This  simple  equation  is  in  good  accord  with  experimental  observations 

2  3  4 

as  reported  in  the  technical  literature  ’  *  for  steady  flow  through  an 
arbitrary  size  and  shape  of  duct  provided  the  pressure  difference  is 
small  conpared  with  the  pressure  at  either  end. 
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The  constant  C  is  termed  the  loss  coefficient  for  the  duct  shape  and 
has  been  measured  for  a  very  wide  range  of  duct  shapes  and  tabulated 
in  the  technical  literature.  Note  that  C  ordinarily  represents  the  loss 
in  total  or.  stagnation  pressure,  measured  in  units  of  the  dynamic  pressure 
between  the  two  points  of  a  duct.  In  most  cases  this  loss  is  not  due  to 
a  frictional  resistance  but  to  an  inertial  resistance  to  lateral  flow  at 
changes  in  area  or  at  bends  in  the  duct. 

There  is  also  experimental  evidence  that  the  loss  coefficient  remains 
constant  for  steady  flow  even  with  large  pressure  differences.  For  large 
pressure  differences,  the  density  cannot  be  assumed  as  constant.  Instead 
we  consider  that  the  flow  is  adiabatic  (neglecting  friction)  and  steady, 
whence  the  kinetic  enthalpy  of  the  air  must  be  constant  along  the  duct,  or 


and,  if  2P  =  p,  +  p  ,  the  flow  equation  becomes 


where  T  is  the  temperature  and  R  is  the  gas  constant. 

This  equation  may  be  compared  with  the  experimental  data  for  the 
steady  flow  of  air  through  a  thin-edged  orifice.  The  experimental 
equation  of  Reference  5  for  subcritical  flow  is  identical  with  Equation 
[16]  if  c  is  taken  as  2.70.  The  experimental  equation  for  supercritical 
flow  is  within  3  percent  of  Equation  [16]  in  the  range  of  pressure  ratios 
up  to  P2/P1  =0.1.  This  range  extends  well  into  the  region  of  supersonic 
flow.  Since  an  equation  of  this  form  is  valid  for  steady  flow  through  a 
thin-edged  orifice,  it  seems  plausible  that  it  should  also  be  valid  for 
steady  flow,  even  with  large  pressure  differences,  through  any  duct  shape, 
because  the  arbitrary  duct  can  be  considered  as  a  series  arrangement  of 
thin-edged  orifices . 
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It  is  convenient  to  define  the  combination  a  =  A  //c  as  the 
"effective  area"  of  the  duct.  It  depends  on  both  the  shape  and  size  of 
the  duct  segment  and,  incidentally,  on  the  direction  of  flow.  The  effect¬ 
ive  area  of  a  complex  duct  system  can  be  calculated  from  the  loss  co¬ 
efficients,  or  effective  areas,  of  its  parts.  Thus,  for  several  duct 
sections  in  series,  where  the  flow  rate  is  the  same  in  all  sections,  the 
effective  area  of  the  combination  is  given  by 


1 


2 


a 


[17] 


While  for  several  sections  in  parallel,  where  the  terminal  pressures 
are  common  to  all  sections,  the  effective  area  of  the  combination  is 

a  =  £  a. 

i  1  [18] 

It  is  precisely  because  the  loss  coefficients,  or  effective  areas, 
can  be  estimated  so  easily,  or  determined  experimentally  if  necessary, 
that  an  attempt  was  made  to  formulate  the  transient  flow  equations  in 
terms  of  these  loss  coefficients. 

An  illustration  of  the  practical  calculation  of  the  effective  area 
and  its  use  in  predicting  the  peak  chamber  pressures  are  given  in  the 
Appendix.  As  a  rough  guide,  it  can  be  expected  that  the  effective  area 
of  a  duct  is  equal  to  half  of  the  area  of  the  smallest  section. 

RANGE  OF  APPLICATION 

Since  the  analysis  is  not  based  on  the  basic  flow  equations  for  a 
gas,  the  range  in  which  the  analysis  is  valid  must  be  determined  by  a 
careful  comparison  of  experimental  results  with  calculated  values.  A 
review  of  some  of  the  assumptions  shows  where  deviations  can  be  expected. 

First,  the  numerical  results  are  valid  only  if  the  external  blast 
pressure  varies  in  the  ideal  way  prescribed  by  Equations  [3]  and  [4]. 

This  variation  is  obviously  a  poor  approximation  in  many  cases  of 
practical  interest.  For  example,  the  blast  pressures  can  be  greatly 
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changed  by  the  presence  of  nearby  reflecting  surfaces  or  diffracting  bodies. 
In  that  case,  it  would  be  necessary  to  modify  Equations  [3]  and  [4]  accord¬ 
ingly,  but  the  resulting  equations  could  still  be  solved  numerically  if  the 
proper  modifications  were  specified.  Equation  [1]  also  assumes  that  the 
incident  blast  wave  has  no  appreciable  velocity  component  into  the  duct. 

2/ 

Otherwise  it  would  be  necessary  to  replace  p^  in  Equation  [1]  by  p^  +  Pv^  '2 
where  v^  is  the  approach  velocity. 

There  are  also  some  obvious  approximations  in  the  boundary  conditions 
for  the  terminal  chamber.  The  isentropic  assumption  neglects  heat  flow 
to  the  walls  and  irreversible  processes  in  the  duct,  but  this  is  estimated 
to  have  a  minor  effect  on  the  peak  chamber  pressure  so  long  as  the  peak 
is  less  than,  say,  2  atm.  It  is  also  assumed  that  the  flow  into  the 
chamber  is  in  immediate  equilibrium  with  the  air  in  the  chamber.  In 
reality,  large  currents  and  large  pressure  gradients  may  exist  in  the 
terminal  chamber,  and  the  peak  pressure  which  is  computed  would  only  be 
a  fair  estimate  if  the  relaxation  time  of  these  currents  were  small  com¬ 
pared  with  the  fill  time  of  the  chamber. 

The  major  approximation  is  probably  due  to  the  lurapimg  assumption. 

This  assumption  immediately  implies  that  we  cannot  deduce  the  pressure 
distribution  in  the  duct  itself.  The  flow  in  the  duct  could  be  unstable 
for  certain  terminal  conditions  whereas  the  analysis  tacitly  assumes 
stability.  The  flow  pattern  in  the  duct,  the  appearance  and  location 
of  shock  fronts,  etc.,  will  all  depend  on  the  instantaneous  terminal 
conditions  and  whether  the  local  flow  is  subsonic  or  supersonic,  but  the 
use  of  a  single  lumped  effective  area  and  single  lumped  longitudinal 
inertia  term  obviously  disregards  these  complexities.  One  justification 
for  this  is  that,  in  practical  cases,  the  major  contribution  to  the 
effective  area  is  due  to  the  entrance  geometry  of  the  duct  and  to  the 
exit  geometry,  and  both  of  these  are  analogous  to  thin -edged  orifices 
for  which  the  effective  area  is  valid  experimentally.^  The  lumped 
inertia  term  is  difficult  to  evaluate  because  it  is  uncertain  how  much 
of  the  total  air  in  the  duct  participates  in  the  average  flow.  Fortunately, 
so  long  as  the  duct  is  small  compared  with  the  chamber  volume,  the  peak 
chamber  pressure  is  rather  insensitive  to  the  precise  value  of  the 
inertia  term. 
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APPENDIX  -  SAMPLE  CALCULATION 


As  an  illustration  in  the  use  of  this  analysis ,  we  consider  the  duct-chamber 
system  shown  in  Figure  4  and  calculate  the  peak  pressure  transmitted  to  the 
terminal  chamber  when  an  ideal  blast  wave  with  a  peak  overpressure  of  20  psi 
and  a  duration  of  0.8  sec  passes  over  the  entrance  to  the  duct. 


We  first  compute  the  effective  area  of  the  duct.  Consider  the 


four  sections,  labeled  (1)  to  (4)  in  Figure  4,  at  which  there  is  a  change 
in  area  or  direction,  and  tabulate  the  loss  coefficient  and  section  area 
for  each.  The  loss  coefficient  is  determined  from  Tables  2  and  3  of  Ref¬ 
erence  2.  We  next  include  the  effective  loss  coefficients  due  to  friction 
and  assume  that  it  is  sufficiently  accurate  to  take  these  as 


_  1  length  of  duct 
f  39  diameter  of  duct 


[19] 


as  is  customarily  done  for  steady  flow  at  low  velocities. 


Circular  Duct 
— 2.0  iq.  ft. 


Circular  Duct 


3.5  tq.ft. 


\ 


Chombtr 
V»760  cuff. 


Figure  4  -  Duct-Chamber  System 
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TABLE  1 


Computed  Values  for  Duct-Chamber  System  of  Figure  4 


4 

Section 

Area 

A 

sq  ft 

Loss 

Coefficient 

c 

C 

A2 

Entrance  (1) 

2.0 

0.85 

0.212 

Bend  (2) 

2.0 

1.30 

0.322 

Expansion  (3) 

3.5 

0.50 

0.041 

Exit  (4) 

3.5 

1.00 

0.082 

Friction  in 
small  duct 

2.0 

0.64 

0.160 

Friction  in 
large  duct 

3.5 

0.12 

0.010 

For  the  entire  duct 


1  ci 

—7  =  E  -4  =  °*827 


A2 

1 


and  the  equivalent  area  a  =  1.10  sq  ft  which  is  55  percent  of  the  area  of 
the  smallest  section.  The  system  fill  time,  in  nondimens ional  units,  is 

0.845V 


T  =  ' 


a  a  e 


=  0.58 


where  a  is  the  initial  sound  speed  in  the  duct. 
The  inertial  ratio  is 


k  = 


0.22 


Now  from  Figure  3  when  the  peak  blast  overpressure  is  1  atm^T  =  0.58, 
and  k  =  0.22,  the  peak  overpressure  in  the  chamber  is  about  61  percent  of 
the  peak  blast  overpressure.  While  for  a  peak  blast  overpressure  of 
2  atm,  t  =0.58,  and  k  =  0.22,  the  peak  overpressure  in  the  chamber  is 
about  54  percent  of  the  peak  blast  overpressure.  Hence,  for  a  peak  blast 
overpressure  of  20  psi  (1.34  atm),  the  peak  chamber  overpressure  would  be 
about  59  percent  of  the  peak  blast  overpressure,  or  12  psig. 
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